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Abstract. Fractional Ornstein-Uhlenbeck process of the second kind 
(fOUa) is solution of the Langevin equation dX t = -6X t df + dY t (1) , 9 > 
with driving noise y t (1) := J* e~ a dB as ; o* = HeTT where B is a 
fractional Brownian motion with Hurst parameter H £ (0, 1). In this 
article, in the case H > ~, we prove that the least squares estimator 
6 T introduced in [0, Statist. Probab. Lett. 80, no. 11-12, 1030-1038], 
provides a consistent estimator. Moreover, using central limit theorem 
for multiple Wiener integrals, we prove asymptotic normality of the 
estimator valid for the whole range H £ (§, 1). 
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1. Introduction 

Assume B = {Bt}t>o is a fractional Brownian motion with Hurst parameter 
H £ (0,1), i.e. a continuous, centered Gaussian process with covariance 
function 

R H (t,s) = l{t* H + s™-\t-s\ 2H }. 

Consider Gaussian process Y t {1) := f*e- s dB as with a t = He^. The frac- 
tional Ornstein-Uhlenbeck process of the second kind X with initial value 
Xq is the solution of the Langevin equation 

(i.i) dx t = -ex t dt + dY t (1 \ x = x , e> o. 

The terminology "of the second kind" is taken from Kaarakka and Salminen 
[8]. The motivation behind the process X is that it is related to Lamperti 
transformation of fractional Brownian motion. For more detailed informa- 
tion on fractional Ornstein-Uhlenbeck processes, see Subsection 12. 11 



An interesting problem in mathematical statistics is to estimate the un- 
known parameter 9 based on continuous observation of the sample paths 
of the process Xt, t £ [0,T]. In the case of fractional Ornstein-Uhlenbeck 
process of the first kind (fOUi ), that is replace driving noise with frac- 
tional Brownian motion B in the Langevin equation (jl.ip . when H > i, Hu 
and Nualart (7j suggested an estimator so called the least squares estimator 

l 
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defined as 

(1.2) 9t = -4-^ 

where the stochastic integral is interpreted as Skorokhod integral (see Sub- 
section [22] for definitions and notations). The motivation comes from the 
following heuristic argument. The least squares estimator is obtain by min- 
imizing the function 

r T 

^ / \X t + 9X t \ 2 dt 



o 

which leads to the solution (|1 .2j) . In the same setup as fractional Ornstein- 
Uhlenbeck process of the first kind, the maximum likelihood estimator of 
drift parameter 9 is found by Kleptsyna and Le Breton [9| and its strong con- 
sistency is proved. Estimation of drift parameter 9 for fractional Ornstein- 
Uhlenbeck process of the first kind, in the non-ergodic case, i.e. 9 < 0, is 
studied by Belfadli, Es-Sebaiy and Ouknine in In this paper, we con- 
sider the least squares estimator 9t in the setup of the fractional Ornstein- 
Uhlenbeck process of the second kind. The key point is the Lemma 13.11 
Using this lemma, we replace the Gaussian noise with an equivalent (in 
distribution) noise and do computations in an equivalent model. 

The paper is organized as follows. In the section [21 we give a detailed infor- 
mation on fractional Ornstein-Uhlenbeck processes and Malliavin calculus 
for factional Brownian motion. Section [3] is devoted to our main results. In 
the appendix section, we provide some auxiliary computation which is used 
in the proof of main results. 



2. Preliminaries 



It is well known that the classical Ornstein-Uhlenbeck process 

{*7 t ( ^ 0) }t>0 with init: 
of the Langevin SDE 



{U t 2 ' ° }t>o with initial value £o can be constructed as the unique solution 



(2.1) dU^' io) = -9U^ o) dt + dW t , t > 0, 

where 9 > and W = {Wt}t>o is standard Brownian motion. The solution 
of SDE (|2.ip can be expressed as 

U &*>) :=e -<*(z 0+ /Vdw,y i>o. 



o 



Let W denote a two sided Brownian motion defined as 



(2.2) W t 



W t for t > 
for t < 0, 



where w[ ^ = {W^ ;t > 0} is another Brownian motion initiated at 
and independent of W. The selection £o = fl^ e es dW s leads to the unique 
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solution given by 

:= f e- fl(( - s) d^, t > 0. 



It is a stationary, continuous Gaussian process with covariance function 
Cw(U$*\u®) = ±e- e \ t -'\, Vs,t>0. 

On the other hand, it is known that the classical Ornstein-Uhlenbeck process 
can be reconstructed from Brownian motion W by Lamperti transfor- 
mation. For a > 0, define the process 

:= e- 9t W ae 2et, ten. 
Then is a stationary, Gaussian process with covariance function 

Since, finite dimensional distributions of a Gaussian process are completely 
characterized by its mean and covariance functions, therefore with a = 

we have Z^ where stands for equality in law. For more 

information on Lamperti transformation and related topics we refer to the 
book [UJ. 

2.1. Fractional Ornstein-Uhlenbeck processes. In this subsection, we 
briefly introduce fractional Ornstein-Uhlenbeck processes. The main refer- 
ences are [5] and [8]. We mainly focus on fractional Ornstein-Uhlenbeck 
process of the second kind that is the core stochastic process of the article. 
To obtain fractional Ornstein-Uhlenbeck processes, replace the Brownian 
motion W with the fractional Brownian motion B in Langevin equation 
(|2.ip . Consider the following stochastic differential equation 

(2.3) duf' io) = -duf' io) dt + dB t 

with initial value £o> The solution of the following SDE can be expressed as 

(2.4) Uf^ = e~ et {i Q + f\ Bs AB t 

Note that the stochastic integral is understood as a path-wise Riemann- 
Stieltjes integral. Let B denote a two sided fractional Brownian motion (see 
[22]). The selection 

£ := / e 0s dB s 



for the initial value £o leads to an unique stationary, Gaussian process 
of the form 



(2.5) 



J — oo 
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Definition 2.1. We call the process given by (|2,4p a fractional 

Ornstein- Uhlenbeck process of the first kind with initial value £o • The pro- 
cess U( H ^ defined in (j'2.5p is called stationary fractional Ornstein- Uhlenbeck 
process of the first kind. 

Remark 2.1. It is shown in [5] that the covariance function of the stationary 
process decays like a power function, so it is ergodic and for H £ (|, 1), 
it exhibits long range dependence. 

Now, we define a new stationary, Gaussian process by means of Lam- 
perti transformation of fractional Brownian motion B: 

X[ a) :=e- at B at , teM, 

where a > and at = —e~H. We aim to represent the process as 
solution of a Langevin type SDE. For this reason, consider the process Y a 
defined via 

Y t {a) := fe- as dB as , t>0. 
Jo 

The stochastic integral is understood in as path-wise Riemann-Stieltjes in- 
tegral. Using the self-similarity property of fractional Brownian motion one 
can see that the process y( Q ) satisfies in the following scaling property 



(2-6) {Y$h>o l = {*- H Y}\>o. 

Using y(°) the process X^ a > can be viewed as the solution of the following 
Langevin type SDE 

dx[ a) = -aX[ a) dt + dY t {a \ 
with random initial value X [ a) = B ao = B H/a ~ JV(0, (f ) 2H ). 

Inspired by the scaling property (j2.6|) . we consider the following Langevin 
equation with Y^ 1 ' as the driving noise: 

(2.7) dX t = -6X t dt + dY t {1 \ 6>0. 
The solution of the SDE (|2.7|) is given by 

(2.8) X t = e~ st (xo + f e Ss dF«) = e~<* (x + f dB a ,\ 

with a = 1 in at- Notice that the stochastic integral is understood as path- 
wise Riemann-Stieltjes integral. The selection Xo = ^~^ s dB as leads 
to an unique stationary, Gaussian process 

(2.9) U t = e- et f e^ s dB as . 

Definition 2.2. We call the process X given by (|2.8|) a fractional Ornstein- 
Uhlenbeck process of the second kind with initial value Xq. The process U 
defined in (|2.9|) is called the stationary fractional Ornstein- Uhlenbeck process 
of the second kind. 
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Proposition 2.1. [8] The covariance function of the stationary process U 
decays exponentially and has short range dependence. 

2.2. Malliavin calculus with respect to fractional Brownian mo- 
tion. In this subsection, we briefly introduce some basic facts on Malliavin 
calculus with respect to fractional Brownian motion. We also recall some 
required results related to fractional Brownian motion which we need for the 
proof of our main theorem. The main references are pQ , [13] and [11] . 

Assume B = {Bt}t£[o,T} is a fractional Brownian motion with Hurst param- 
eter H S (0, 1) and covariance function Rn(t,s). It is well-known that the 
covariance function Rh can be represented as 

ptAs 

R H (t,s)= / K H (t,u) K H (s,u)du 
Jo 

for a Volterra-type square integrable kernel Kh- In the case H > i, the 
kernel Kh has a simple expression given by 

1 ft 3 1 

K H (t,s) = c H s2~ H (u - s) H -2u H ~2du 

J s 



n/ 2i?r(f-if) \5 



where ch = (H — o)(r=— — nTT . I . Consider the set £ of all step 

functions on [0, T\. We define the Hilbert space T~L as the closure of E with 
respect to inner product 

( 1 [o,t], 1 [0,8])H = RH(t,s). 
Then, the mapping lr ^ i->- Bt can be extended to an isometry between 
Hilbert space T~L and Gaussian space %\ associated with fractional Brownian 
motion B. We denote this isometry by y i — > B((p). It is known that when 
H = 2, we have % = L 2 [0, T], whereas for H > ^, the elements of % may be 
not functions but distributions of negative order (see [E]). An interesting 
subspace \H\ containing only functions is the set of all measurable functions 
<p on [0, T] such that 

\\<p\\f H \ : = a H J J \ip(s)Mt)\\t - s\ 2H ~ 2 dsdt < oo, 

where an = H(2H — 1). Notice that for any two measurable functions 
(p,ip & \~H\, we have 

(2.10) m{B{tp)B{ij))) = {<p, il)) u = a H f C <p(s)il>(t)\t - s\ 2H ~ 2 dsdt. 



Notice that when H > \, we have the inclusions L 2 [0, T] cLh[0, T] C \7i\C 
H. Consider the linear operator K* H between £ and L 2 [0,T] defined by 

(K H ^)(s) = f\(t)^{t,s)dt. 

Notice that the operator K H is an isometry between £ and L 2 [0,T] that 
can be extended to the Hilbert space T~L. Moreover, the processes W = 
{ w t}te[o,T] given by 

(2.H) W t :=B{{K* H r\l m )) 
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defines a Brownian motion. The fractional Brownian motion B and Brow- 
nian motion W relate through the integral representation 

(2.12) Bt= f K H (t,s)dW s . 

Jo 

Consider the space S of all smooth random variables of the form 

(2.13) F = f(B(ipi), ■ ■ ■ , B((p n )), ipi, - • • ,ip n £ H, 

where / G C^°(IR n ). For any smooth random variable F of the form [2~T3l 
we define its Malliavin derivative = D as an element of L 2 (tt;n) by 

n 

DF = Y i a i f(B(ip 1 ) ) ---,B(ip n ))<pi. 

i=l 

In particular, DB t = l [0)t] . We denote by B^ = B 1 ' 2 the Hilbert space of 
all square integrable Malliavin derivative random variables as the closure of 
the set S of smooth random variables with respect to norm 

\\F\\l 2 =m\F\ 2 +m{\\DF\\D. 

The transfer principle connects the Malliavin operators of the both processes 
B and W . This is the message of the next proposition. 

Proposition 2.2. p] For any F G B^ 2 = B 1 ' 2 , 

K* H DF = D {W) F, 

where D^ w ^ denotes the Malliavin derivative operator with respect to Brow- 

1 2 

nian motion W , and the corresponding Hilbert space. 

The divergence operator denoted by 5b = S is the adjoint operator of Malli- 
avin derivative operator D. For any %-valued random variable u £ L 2 (£2; H) 
belongs to the domain Dom<5 of the divergence operator 5, then the random 
variable 5(u) is defined by duality relationship 

JE(FS(u)) = JE(DF, u) H , VF e B 1 ' 2 . 

Also, an element u € L 2 {Q,;T-L) belongs to the domain Dom<5 if and only if 

\TE(DF,u) H \ <c u \\F\\ L 2 

for any F G B 1,2 , where c u is just a constant depending on u. The divergence 
operator 5 is also called Skorokhod integral. It is known that in the case 
of Brownian motion it coincides with Ito integral for adapted integrands. 
Hereafter for any u G Dom<5, we denote S(u) = Jq u t 5B t . 

Following pQ one can develop a Malliavin calculus for any continuous Gauss- 
ian process G of the form 

G t = I K(t,s)dW s 
Jo 

where W is a Brownian motion and the kernel K, < s < t < T satisfying 
su Pte[o,T] Jo K(t, s) 2 ds < oo. Consider the linear operator K* from £ to 
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V 



L 2 [0,T] defined by 

(K*<p)(s) = <p(s)K(T, s) + [ T [<p(t) - <p(s)] K(dt, s). 

J s 



The Hilbert space T~L generated by covariance function of the Gaussian 
process G can be represented as U = (A'*)~ 1 (L 2 [0, T]) and B^ 2 (%) = 
(K*)- 1 ($#{I?lp,T\)). For any n > 1, let Jf n be the nth Wiener chaos 
of G, i.e. the closed linear subspace of L 2 (f2) generated by the random 
variables {H n (G((p)) , (p G H, ||(^||% = 1}, and H n is the nth Hermite poly- 
nomial. It is well known that the mapping I^{tp® n ) = n\ H n (G((p)) provides 
a linear isometry between the symmetric tensor product 7-L® n and subspace 

Definition 2.3. [1] We say that the kernel K is regular if for all s G 
[0, T), K(-,s) has bounded variation on the interval (s,T], and 

r T 

/ \K\((s,T],s) 2 ds < oo. 
Jo 

For regular kernel K, put K(s + , s) := K (T, s) — K((s, T], s). For any ip G £, 
define the seminorm 

ll^llL=^ T ^) 2 ^ + ,s) 2 d S + ^ T ^ T |^)||A|(dt, S )) ds. 

Denote by T-Lkt the completion of £ with respect to seminorm || ,\\kt- The 
following proposition establishes the relationship between path-wise integral 
and Skorokhod integral. 

Proposition 2.3. PQ Assume K is a regular kernel with K(s + ,s) = and 
u is a process in (Hkt)- Then the process u is Stratonovich integrable 
with respect to G and 

[ u t dG t = [ u t SG t + f ( [ D s u t K(dt, s)) ds. 

JO Jo Jo \Js J 

The next proposition is taken from [T2] and provides a central limit theo- 
rem for a sequence of multiple Wiener integrals. Let M(0,a 2 ) denote the 
Gaussian distribution with zero mean and variance a 2 . 

Proposition 2.4. Let {F n } n >i be a sequence of random variables in the q- 
th Wiener chaos, q>2, such that lim„_> 00 ~\E(F 2 ) = a 2 . Then the following 
statements are equivalent: 

(i) : F n converges in distribution to M(0,a 2 ) as n tends to infinity. 

(ii) : HDFjill 2 ^ converges in L 2 (f2) to qa 2 as n tends to infinity. 

The next proposition will be used to bridge our original Langevin equation 
(|2.7p to a new SDE that plays the key role in our computations. 

Proposition 2.5. [3j Let B = {Bt}te[o,T] be a fractional Brownian motion 
with Hurst parameter H G (^,1). Suppose f : [0,T] — > IR be a strictly 
positive, absolutely continuous and f is locally square integrable function. 
Then there exists a locally square integrable kernel L such that 

B rtmi*. '=[ !{8)AB ' + t ( [ aJ t(u,s)L {a ,r )A ,yuyw r 
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where Brownian motion W is given bv \2.11[ Moreover the kernel L satisfies 
in 

cf>(t, s)=£ f{t) (l\ 2 {t - r) H ~i L(r, s)dr 

and 

*(*,») = ^ff ( /'/(»)**'./'/(<>)**')/(«)*/(»)* 



3. Main Results 



For the rest of the paper, we assume that H > |. Let X = {-X~t}tg[o,Tl be 
the solution of the Langevin equation (|2.7p . Assume Xo = and 9 > 1 (this 
is a technical assumption which we need in the proof of the main result). 
Then the solution can be expressed as 

(3.1) X t = fe-^dYP. 

Jo 

The next theorem gives the main result of the paper. 

Theorem 3.1. The least squares estimator 9t given by (|1.2|) is weakly con- 
sistent, i.e. 

9 T ^9 

in probability, as T tends to infinity. 

The next theorem provides the rate of convergence for the least squares es- 
timator towards the Gaussian distribution. Let B(x,y) denote the complete 

Beta function and '-^T stands for convergence in distribution. 

Theorem 3.2. For the least squares estimator 9t given by (|1.2|) . we have 

Vt (e T - e\ Af(0, a 2 (9, H)) 
as T tends to infinity, where 

J2l 



a\9,H) 



H 2 B((9- 1)H + 1,2H- I) 



V 



e h — e h\\\ 



X 2 / e -6y e -e\z-x\ e (l-^)(x+y+z) 
J[0,oo) 3 

Remark 3.1. It is worth to mention that the central limit theorem \3.2\ holds 
for whole range H £ (^,1), whereas for the fractional Ornstein-Uhlenbeck 
process of the first kind we have restriction H G [|, |), see [7]- 

To prove Theorem 13.11 we need a series of lemmas. We start with the 
following lemma that provides an alternative stochastic differential equation 
of the Langevin equation ([27 



2H-2 

dxdydz. 



Lemma 3.1. Let Bt = Bt+H — Bh be the shifted fractional Brownian mo- 
tion. Then there exits a regular Volterra-type kernel L, in sense of Definition 
\2.3l so that for the solution of the following stochastic differential equation 
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(3.2) dX t = -6X t dt + dG t , X = 0, 

we have, {^t}te[o,T] {Xt}te[o,T} where the Gaussian process 

G t = J t (K H (t,8)+L(t,sj) dW s 
and Brownian motion W is defined as \2.11\ 

Proof. It is easy to check that the solution of the SDE (|3.2p is given by 

X t = I e-^dGs, te[0,T}. 
Jo 

Moreover, E(X t ) = E(X t ) = 0, and 

e -0(t-u) d y(i) f S e -e(s-v) dY (i) 



t 

I e -e(t-n)-9(s-v) CoVy(1) ( dU) dv y 
JO Jo 
Similarly 

Cov^(t, s) = E (J e - 9{t ' u) dG s jT e - 9{s - v) dG^j 



t r s 

e -9(t-u)-e(s-v) Cov & (du,dv). 
JO 



Therefore, it is enough to show that for the Gaussian processes Y« and G, 
we have 

Cov y( i)(M) = nY s (1) Y t {1) ) = E(G s G t ) = Covg(a,t). 
To see this, consider the function f(t) = e t and apply proposition ^. 51 Then 

(3.3) B at = Bjt ^ dr + B H ^ e r dB r + g(t, r) dW r + B H 

where the Volterra-type Kernel g is given by 

g(t, r) = f e u ( f ^ (u, s)L(s, r) ds\ du. 



ds 

Now, using integration by parts formula, we obtain 
(3-4) 

Y t W = I e- s dB as 1 ^ f e~ s d( f e r dB r + f g(s, r) dW r + B H ) 
Jo Jo \Jo Jo / 

= Bt + J o e^dU* g(s,r)dW r \ 

= B t + e-* J g(t, r) dW r + f* QT g(s, r) dW^j e~ s ds 
Next, using Fubini theorem for Brownian motion W, we see that 

A = j (j g(s,r)dW^J e~ s ds = J j e~ s g{s,r) ds dW r 
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Put 

h{u,r) = J ^ H (u, v)L(v, r) dv, 
and using Fubini theorem for real-valued functions, we obtain 



I e s g{s,r)ds = f e 3 f 

J r J r J r 



e u h(u, r) duds 



r J r 

e u h(u,r) (J e- s ds^j du 

t rt 

t 



h(u,r)du — e \ e u h(u,r)du. 
Now, plug this into A, we infer that 



A= [ [ e~ s g{s,r)dsdW r 
Jo Jo 

= Jo (/ h ^ r ^ du ~ e ~ t J e u h(u,r)du^j dW r :=A 1 -A 2 

where 

(3.5) Ai = J 1^ h(u,r)dudW r , 
and the term A 2 can be written as 

(3.6) A 2 = [ e~* / e u h(u,r)dudW r = [ e~*^(t, r) dW r . 

Jo Jr Jo 

Finally, plug in ([53]) ■ (1331) and (l2J2]l into (HQ}, we obtain 

y« j* (K H (t, s) + J* h(u, r) dv\ dW r . 
Therefore, it is enough to introduce the kernel L as 

(3.7) L(t,r) = J h(u,r)du. 

and then the result follows. □ 

The next lemma gives an alternative form for the least squares estimator in 
terms of the process X and Gaussian noise G. 

Lemma 3.2. For the least squares estimator 9t we have 

(3.8) e T l =9- £ Xt5Gt 



f Q T X?dt' 

Proof. Note that the least squares estimator Qt is actually equal to 



0t = 



10 

On the other hand, we have 



fixfdt 
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X t 5Y t {1) = f T f e~ e ^ SYP 5Y t W = if 1 '(/) 



o J o 

and the function / is given by f(t,s) = i e — *' . Similarly, 

/ T X f SG t = [ T f e~ 9 ^ 5G S SG t = 4(f). 
Jo Jo Jo 

This immediately implies that 

T Xt6Y W law [ X t SG t . 



o 



Moreover, for every i G [0,T] there exists a function g(-,t) so that we can 
write 

X t = lY W (g(-,t)). 
Using multiplication formula for multiple stochastic integrals, we obtain 

X? = \\ 9 (;t)\\ 2 H + ir i) (g(;t)®g(;t)) 

where T~L stands for Hilbert space generated by the covariance function of 
the Gaussian process . Therefore, using Fubini theorem 

J X ? dt = J Q Il5(-,i)ll«dt + / 2 y(1) (J {g(;t)®g(;t)) dtj 

Similarly, in the same lines one can obtain that 

J X? dt = J Q \\9(;t)\\ldt + li (J (g(;t)®g(;t)) dt 

where % stands for Hilbert space generated by the covariance function of 
the Gaussian process G. Note that i)||-H= HffO)^)!!^- Consequently 

T xfdt 1 ^ f Xldt. 



o 



Finally, by Lemma 14.21 in Appendix, we can deduce that for the following 
two dimensional random vectors we have 

\M'\[xUt) la = (]*X t 5Gu[x 2 t ^ ■ 
We conclude that 

Jo Xt SY t j Q X t 5Gt 



tfxfdt tfxfdt 



□ 



In the next lemma, we give an alternative expression of the estimator Ot- 
Lemma 3.3. The least squares estimator 9t can be written as 



Xl fo II ^Xt (K H (dt, s) + L(dt, a)) ds 



^ law \X-t . J° J s 

T ~ " 'f?X?dt tfxfdt ' 

where D^ w ^ denotes the Malliavin derivative operator with respect to Brow- 
nian motion W . 
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Proof. Using Lemma [3.21 Proposition 12.31 and SDE (|3.2[) . we can infer that 
a law„ Jq X t SG t 

f Q T X t dG t - f Q T jf Df ] X t (K H (dt, s) + L(di, s)) ds 

= 9 TfT^Z 

JoXfdt _ 

/ T X t (dXt + 0X t dt) - J T jf £>f°X t (^(dt, s) + L(dt, s)) ds 



Jo 



£x*dt £x?dt 

□ 

In order to prove the theorem, we need the following lemma. The next 
lemma can be used for an alternative estimator of the parameter 9. 



Lemma 3.4. For the process X given by (|3.ip and consequently for the 
process X , as T tends to infinity, we have 

I r V 2 A < x m - i)h** 

T .„ 

almost surely and in L? . 
Proof. Note that for every t > 0, we have 

(3.9) U t ~ 



f T X 2 t dt -» {2H l )R2H B ({9 - 1)H + 1,2H- 1) 
Jo 9 



f e^ s dB as =X t + e~ s % 

J — oo 



where £ = J*^^ ^ s dB as . For the stationary Gaussian process U we have 
(see Proposition 3.11 of [8]) 

E (U t U ) = 0(e- et ) -> 0, 

as t tends to infinity. Hence, U is ergodic. Therefore, the ergodic theorem 
implies that 



^£ufdt^E(U 2 ) 



10 

as T tends to infinity, almost surely and in L 2 . It is straightforward from 
(pT5j) to check that 

(3.10) i^ T X t 2 dt^E([/ 2 ). 

On the other hand, using change of variable we see that 

r-t 



U t = e~ et [ e^dB as 

J — oo 

= H -(e-i)H e -<n f at s (e-i)H dB 
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Therefore, using formula (|2.10p . we obtain 
(3.11) 

E([/ ) 2 = H-W-V H E QH s^ H dB a \ 



H -W-W H(2H _ X) / / s (ff-l)H t (ff-l)H| s _ fl2ff-2 dsdt 



JO 

= ( 2H -^ H2H B {{9 _ 1)H + i, 2if - 1) 
Substituting (|3.1ip into (|3. 10)) yields the result. □ 

Proof of Theorem I3.il From lemma (|3.3|) . we can write 9t = h + h- By 
Lemma 13.41 we know that ^ Xf dt converges almost surely and in I? , as 
T tends to infinity to a constant. Then, it suffices to study the almost sure 
convergence of the numerators of I\ and Ii- 

Step 1: I\ — > almost surely as T tends to infinity. From Lemma (|4.ip in 
Appendix and the assumption 6 > 1, we deduce that almost surely 

X 2 

(3.12) lim — f = 0. 

Step 2: I2 — )• almost surely as T tends to infinity. 

First note that the Malliavin derivative Ds^ Xt with respect to Brownian 
motion W is given by (see Section 2.2) 

Df )X t = J* e*™ f^(u, s) + h(u, s)\ du. 



Also, 



K H {dt, s) + L(dt, s) = (^-(t, s) + h(t> s )) d * 

= (^ L (t,s) + e- t <f>(t,8)\ dt 

_ f dK H //■*«, f s u , 

= e — — / es du, / e» du e^ean dt. 



at 

Denote the numerator of the term I2 with i^ p . Therefore, 



IT = [ D^Xt (K H (dt, s) + L(dt, s)) ds 



= — / / e { >e — — ( / en du, ea du I e» ean 

x e — — / e» du, / e« du eae^a df ds. 

Now, by Fubini theorem for real-valued functions and the Lemma (|4.3p in 
Appendix, we obtain 
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1 2 



1 



T i-t 



J 



T 

tAu QK 



11 







Off 

T 



ft 

T ft 



JO 



e» dv, J e~H dv 
«M e -« e f e -t e l 



o 



air 



du 
&h dv 



o 





2_ff-2 



en dv, / e» dw e^dsdudi 



du dt 



We make a change of variables p = ea dv and (/ = J* * ea dv. Hence, 



J 2 



T 



aj-~H fQ 



1 + i 



-#(1+0) 



1 + ¥ 



p \J5r(e-i) 



|p - g 



2i/-2 



dpdg 



^(fl-i) 



p-q 



\2H-2 



dpdg 



H 2H pot l rl 



t'u 



-1)H, 



l-y) 2H - 2 dydq 



As a result, 



I% p — > (2# - l)ff 2H 5 ((0 - 1)# + 1, 2# - 1) 
and the claim follows. 



□ 



Proof of Theorem \3.2l Using Lemma 13.21 we have 



VT(6 T -e 



law 



Vti? 



2 l2< 



'i a -e\t-s 



Ft 



fixfdt 



1 r± Y 2 Hf 
;o u " t Jo ^* ut 

where i*r stands for the double stochastic integral 



1 



F T = — 7 2 G ( 



Therefore, taking into account Lemma 13.41 it is enough to show that the 
sequence {Ft} for T = 1, 2, • • • converges in law to a Gaussian distribution. 
To show this, we use Proposition 12.41 We have 

.,2 r . 

-9(\u 1 ~v 1 \ + \u 2 ~v 2 \) e (jj-l)(u 1 +v 1 +u 2 +v 2 ) 



"2 



[0,7? 

e»dx 









2H-2 


rv 2 


/ x 

/ eHdx 






Jo 




Jo 



e^dx 





2H-2i 


/ e^dx 




Jo 





duidu2d-yidf2 



_ a H tt4H-4t 
■ 2T T, 



where 



I T 



[0,T] 4 



(|ui— «l| + |tt2— Val) „(-^-l)(Ml+fl+M2+f2) 



duidu2dt>idt>2. 





"2 "1 


2H-2 


v 2 


2H-2-1 


X 


e h — e h 
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Taking the derivative with respect to T, change of variables x = T — u±,y = 
T — U2, z = T — v±, and taking limit as T tends to infinity, we obtain that 



Jim ^ = 4 / e -<fo e -»l*-*l e (i-*)(*-H<+*) 



T->oo dT 



[0,oo) 3 

2.H-2 



X 



e -ff — e -tf 1 — e «1 



dxdydz 



and this integral is finite. Indeed using obvious bound e 9 \ z x ' < 1 and 
change of variables = a, e~ «" = 6, e~"& = c, we infer that, it is smaller 
than 



C{H)B{1 — H, 2H — 1) / (afc)-^^ - a| 2 "- 2 dad& 

J [0,1] 2 

= C(H)£(1-H,2#-1)[ /* / + f 1 (ab)~ H b eH \b - a\ 2H ~ 2 dadb 

Ljo jo jo A 
= C{H)B{1 -H,2H- 1) ( J x + J 2 ) . 

For example, for the term J 2 we have 

J 2 = b m - l db C w^ H (l - w) 2H - 2 dw 
Jo Jo 

= B(l + (6- 1)H, 2H - 1) f a eH ~ l da < 00. 

Jo 

With similar computation (see [2], appendix), one can show that 

2 



I^FtI^-ieiia^tII! 



as T — > 00. 



IE 

Taking into account that 

lim E||L> s F r || 2 .= 2 lim JE(f£), 

we complete the proof. □ 

Remark 3.2. Note that if one can replace Skorokhod integral with path-wise 
Riemann-Stieltjes integral in the formula of least squares estimator, then we 
can obtain the new estimator 

-r fn x t dx t 
tfxfdt 

almost surely by Lemma \3.4\ and \4-l\ in Appendix, as T — )■ 00 . 

4. APPENDIX 

Lemma 4.1. For the stationary Gaussian process U given by (|3.9p and any 

a > 0, as T tends to infinity, we have 

U T 

— > almost surely. 

Proof. Using Proposition 3.11 of [8], there exists a constant C(H,8) such 
that 

rat f-ao 

(xy)^ H \x-y\ 2H - 2 dxdy 



(rat r° 
Jo Jo 



16 



AZMOODEH AND MORLANES 



and the constant C(H,9) = H(2H - l)H 2H ^~ e \ Denote the term inside 
parentheses by I(t). Then with some direct computations one can see 



a t ' 



I(t)= / at (xyf-^(y-xr-'dxdy 
Jo Jo 

+ [ X (xyf-^ H (x-y) 2H - 2 dydx 
Jo Jo 

+ r f {xyf^ H {y-xf H - 2 dxdy 
Jo J^-y 

= ^B{(6 - 1)H + 1, 2H - 1) + ^(af H - af H ) / 1 z^(l - zf^dz. 
Therefore, as t — > + , we have 



(4.1) 



c(t) = (2ff l)H 2 _ i)R + l2R _ 1} 

u 

- (2H - l)H 2H xtx I'' z^ H (l - z) 2H ~ 2 dz + o(t 



Let Bi nc (p, q, x) and 2Fi(a, b, c; x) stand for incomplete Beta function and 
Gauss Hyper geometric function respectively (see [17] for definition). Using 
a relation between these two functions ([IT], page 289), we obtain that as 
t — > + , we have 

(4.2) 

t z^ H (l - z) 2H - 2 dz = B iQC (2H -1,(9- 1)H, 1 - 



JTTT-2H 

e^\a t - a ) 2H ^ ^{{9 + l)H, 1, 2H; 1 - ^) 

t 2H-l +o{t 2H-ly 



2# - 1 

Substituting (|4.2p in (|4.1|) . we obtain that as i — >• + , we have 

= (2g -l)H 2H B ^ e _ +li 2H-l)- Ht 2H + (t 2 "). 
9 

Now the claim follows from Theorem 3.1 of [15]. □ 

Lemma 4.2. Sei C7 = I% m (/), F = lJf W ( J Q T (g(-,t)®g(; t)) dtj and U = 
li(f), V = 1$ (/ T (</(-, *)®<7(-,t)) dt). Tften 



E(C/ mi y m2 ) = E (^[/ mi y m2 J Vmi,m 2 > 1. 
Proof. It immediately follows from Corollary 7.3.1 of [T3]. □ 
The next lemma is used in study of the term I% p . 
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Lemma 4.3. 



sAt 



dK 



H 



at 



e H dv, / es dv 



I e» du 

2H-2 







#(2F - 1) 



o 



n 



Proof. Put Vt = B^^^. Note that we can also represent this random 
variable as Wiener integral in the following way: 



K H I e~H dv, \ e~H dv) e^n dW u 



Vt 



where Brownian motion W is given by (|2.1ip . Therefore 

(4.3) 

E (V 8 V t ) 



sAt 



K H 

On the other hand, 



en dv, / en dv) K H I e» dv, / e«dw e» du 



(4.4) 



E(V s V t )=E(B rs _v rj _B l . t 



Rh 



So elT dv fo elT dv 

e^ dv, I en dv) . 



o 



With differentiating respect to s and t of the right hand side of (|4.4p and 
(|4.3p . we get the desired identity. 

□ 
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